fig = figure;
errorbar(trend,U50(:,3),U100(:,3)-U50(:,2),U50(:,4)-U50(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Gray'));
hold on;
  errorbar(trendp,M50(:,3),M100(:,3)-M50(:,2),M50(:,4)-M50(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Black'));
hold off;
xlim([0 114]);
%ylim([-0.02 0.09]);
orient landscape;
print(fig, [figdir, 'Figure_6_a_50yr' mod_name], '-dpdf','-bestfit');
close(gcf);

fig = figure;
errorbar(trend,U100(:,3),U100(:,3)-U100(:,2),U100(:,4)-U100(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Gray'));
hold on;
  errorbar(trendp,M100(:,3),M100(:,3)-M100(:,2),M100(:,4)-M100(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Black'));
hold off;
xlim([0 114]);
%ylim([-0.02 0.09]);
orient landscape;
print(fig, [figdir, 'Figure_6_a_100yr' mod_name], '-dpdf','-bestfit');
close(gcf);